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Abstract 

The program package XLOOPS calculates massive one- and two-loop Feynman 
diagrams. It consists of five parts: 

• a graphical user interface 

• routines for generating diagrams from particle input 

• procedures for calculating one-loop integrals both analytically and numerically 

• routines for massive two-loop integrals 

• programs for numerical integration of two- loop diagrams. 

The package relies on the application of parallel space techniques. The treatment of 
tensor structure and the separation of UV and IR divergences in analytic expressions is 
described in this scheme. All analytic calculations are performed with MAPLE. Two- 
loop examples taken from Standard Model calculations are presented. The method 
has recently been extended to all two-loop vertex topologies, including the crossed 
topology, graphs with divergent subloops and IR divergent diagrams. This will be 
included in the XLOOPS package in the near future. 
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XLOOPS 



an introduction 



to parallel 



space techniques ^ 



Dirk Kreimei]^ 



The package XLOOPS presented in this workshop rehes on the apphcation of parahel 
space techniques. We introduce these techniques covering the following topics: 

• The generation of integral representations for massive two-loop diagrams. 

• The treatment of tensor structures. 

• The handling of the 7-algebra in this scheme. 

• The separation of UV and IR divergences in analytic expressions. 
We present two-loop examples taken from Standard Model calculations. 

1.1 Introduction 

In last year's conference in Pisa, AIHENP95, David Broadhurst and myself presented 
for the first time results concerning a fascinating connection between knot theory and 
Feynman diagrams.^ Meanwhile this has lead to numerous results in field theory, as well 
as in knot theory and number theory.'^ We are still working in that area, and the results 
are even more intriguing, with the connection between transcendental numbers arising in 
counterterms and knots associated with the Feynman graphs becoming more and more 
specific. David Broadhurst and I just finished a joint two month stay at the University of 
Tasmania, where Bob Delbourgo invited us to continue our collaboration on the subject. 
As a consequence, David had to go back to teach at summer schools directly after our stay 
in Tasmania, so that for him it is not possible to attend this conference. 

Today, I would like to focus on the techniques used in the package XLOOPS presented 
at this workshop. It is based on a very simple idea: that it might be useful to combine 

^Work supported by grant CHRX-CT94-0579, from HUCAM. 
^email: kreimer@dipmza.physik.uni-mainz.de 
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dimensional regularization with the fact that each Green function offers a distinguished 
subspace of spacetime - the space spanned by all its external momenta. This very fact was 
used as a convenient means to define dimensional regularization, for example in Collins 
text book.^ 

Amazingly, a systematic exploration of this idea allows to make progress for the calcu- 
lation of general massive one- and two-loop Green function. The package, as it is presented 
in this conference, focuses on two- and three-point Green functions.^ While Lars Briicher 
and Johannes Franzkowski will report on the package XZO OPS' itself, Alexander Frink will 
report on progress for the scalar two-loop three-point functions, covering planar and non- 
planar topologies, results which were already used in recent Standard Model calculations.^ 



1.2 Parallel and Orthogonal Spaces 

In our approach, we use the fact that any non-trivial Green function furnishes a set of 
external momenta, on which it depends. A n-point Green functions provides in fact n — 1 
independent external momenta pi. This allows to write any loop momentum ki as a sum of 
two covariant vectors, fc^ = fc^y + k'^j_, where fejji has components in a space which is the 
linear span of the external momenta pi, the parallel space, while ki^± is in its orthogonal 
complement, ±_PiiJ. = 0- Dimensional continuation happens in the orthogonal space. 



1.2.1 Two-point functions 



To consider the first specific example, we turn to a massive one-loop two-point function, 
providing a scalar integral of the form 
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:, k±-q = 0. 



(1) 



(2) 
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We straightforwardly integrate the orthogonal space to find expressions which involve 

[^iUx=o]'^, (4) 
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which determine the cut structure of the result. Note that we avoid the Wick rotation 
and Feynman parametrizations altogether. The remaining integration with respect to c?A;|| 
can be done easily. For example, we can reinterprete it as an integral representation of a 
"/^.-function.^ Upon expanding in D — 4, we recover the standard results of perturbation 
theory. 

If we go to the two-loop level, we have five non-trivial integrations. We had two for the 
one-loop case, and for the two-loop case we have one more than twice this number, as there 
is an extra angular integration in orthogonal space, due to the presence of a non-trivial 
scalar product l± ■ k±. It is known that straightforward application of the parallel space 
method gives a two-fold integral representation.^ This integral representation is the basis 
for the two-loop two-point package included in XLOOPS. 

1.2.2 Three-point functions 

Now we have a two-dimensional parallel space. Accordingly, in the one-loop case, the scalar 
one-loop integral has the following form in parallel and orthogonal space variables 

-D-2 , n 

Pi = - kj - kl - ml, 

P2 = {ko + qiflf - (ki + qi^if - kl - ml 

Ps = {ko + q2,o? - {ki + q2,if - kl - ml (5) 

where 

/t/|l — ^0^1 ^i^^j — ~\~ k I ^ 
ei ■ k± = €2 ■ k± = 0, span{ei,e2} = span{qi,q2}. (6) 

There are two non-trivial integrations for the parallel space, assumed to be spanned by two 
external momenta qi,q2, and one for the modulus of k±. One can proceed by explicitly 
using the signature of spacctime in the parallel space: a shift /cq ^ + ki renders all 
propagators linear in the variable ki. The orthogonal space integration confronts us with 
a non-trivial cut structure: the integrand obtains a factor [Pi \i^_^^q]^~ . Due to the shift 
in ko, the cut is confined to a half plane in complex ki space, so that we can use the residue 
theorem for one of the parallel space integrations. The resulting expression can be identified 
as an integral representation of a i?-function or, equivalently, of a hypergeometric.^ 

In the two-loop case, some of the features remain. We are now confronted with four 
integrations for the parallel space, and three for the orthogonal space. As in the two- 
point case, we do the angular integration in the orthogonal space, with respect to z = 
jl^i^, first. This results in a non-trivial cut structure, and the remaining integrations are 
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determined by the demand that as many as possible can be done as residue integrals by 
closing the contour at infinity. A residue contributes only if it is located in the interior of the 
contour, and this results in constraints for the remaining variables. It turns out that these 
constraints determine the remaining domains of integrations to be finite triangles in the 
space of two remaining parallel space variables. At this level we are talking about a four-fold 
integral representation. Two further integrations in the orthogonal space variables 
can be achieved for the planar topology, resulting in a two-fold integral representation, 
by using Euler transformations.® Alexander Frink achieved a similar result for the non- 
planar topologies, which was already used and tested recently.^ He will report on these 
achievements in his talk in detail. 



1.2.3 The general case 

One could extend the method to four-point Green functions (and higher). Wc have not 
included this at this stage, but do not see any conceptual problems to do so if necessary. 



1.3 Tensor Structures and 7- Algebra 

One of the most interesting features of our approach is its treatment of the spin structure of 
Green functions. We altogether avoid the calculation of single tensor integrals (though this 
is possible in XLOOPS, to make contact with more conventional methods), but directly 
calculate the characteristic polynomials of a graph.^ 



1.3.1 Two-point functions 

Again let us go back to a one-loop two-point function to understand the idea. Any Green 
function will deliver polynomial expressions in kg, k'j_ for its numerator. They might come 
distributed over various form factors.^ To handle this situation, we also separate the Clif- 
ford algebra into two orthogonal sub-algebras, according to the splitting into parallel and 
orthogonal spaces. Using basic properties of Clifford algebras, it is then easy to determine 
the characteristic polynomials. A trivial one-loop example is provided by the self energy 
of a fermion, dressed with a massless vector boson in the Feynman gauge 

/■ „ „ ln[k\\J\\ - 7± ■ k± + m]-f^' f [(2 - D)fc||7|| + ml] 

J '^''11 "^"^ p + g)2][P_m2] ' J ^''W ^^^^JkTWW^^' ^ ' 

We have two formfactors, one proportional to 7||, the other one proportional to 1 in spin 
space. In parallel space, the only element of the 7-algebra is = q'^q^'Ju/q'^ = 7-6, while 
the Clifford algebra in orthogonal space is spanned by elements 7^, with q ■ J± = and 
{711 )7±} = 0.^ In the above example we found two very simple characteristic polynomials 
for the two formfactors, (2 — D)k^^ and m. 
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Note that any polynomial expression in kp k^ can be calculated by either reducing 
the expression to some basic scalar integrals, or to trivial massive tadpoles, by using the 
following rules 

P2-P1 + [ml -rnj- q^] 
k\ = k\-Pi- ml, 

Johannes Franzkowski will have further comments on the implementation of these ideas 
at the two loop level in his talk. 



1.3.2 Three-point functions 

Here we have the possibility to reduce polynomial expressions in the variables ko, ki, k^ to 
either two-point functions, or to basic scalar three-point functions. Again, we achieve this 
by solving for the above variables in terms of propagators, masses and external momenta. 
These reductions are already implemented in XLOOPS in the one-loop case, and will be 
available at the two-loop case in the near future. 



1.4 UV and IR Divergences 

We usually do a termwise determination of the UV degree of divergence for our character- 
istic polynomials. We subtract massless integrals to achieve finite integral representations, 
carefully avoiding oversubtractions and thus avoiding the generation of spurious infrared 
singularities.^ A main feature of our approach is that the subtraction is done termwise in 
the characteristic polynomial. This avoids oversubtractions for the non-leading terms in 
the characteristic polynomials, while only the powercounting for the leading terms in the 
characteristic polynomial agrees with the powercounting for the graph itself. 

At the moment, we gain first experience for the subtraction of IR singularities from 
our integral representations in collaboration with Jochem Fleischer in Bielefeld. 

At the two-loop level, we calculate the UV divergences of a Feynman graph analytically, 
and generate well-defined code for numerical integrations.^ 

1.5 Conclusions 

Let me list our main results as follows: 
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• XLOOPS calculates arbitrary graphs in the Standard Model for one- loop two- and 
three-point functions. After specifying a consistent particle content for a chosen 
topology, the user can choose to get an analytic or numeric result. 

• At the two- loop level, XLOOPS' returns an analytic result for the UV-divergent part 
of any two-point Standard Model graph, and automatically generates code for the 
remaining finite part of the whole graph. 

• We have integral representations for all scalar two-loop three-point functions and 
will implement the two-loop three-point case in XLOOPS in the future. 
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2 



XL OOPS — a package calculating 
one- and two-loop diagrams 

Lars Briichei]^ 



A program package for calculating massive one- and two-loop diagrams is introduced. 
It consists of five parts: 

• a graphical user interface 

• routines for generating diagrams from particle input 

• procedures for calculating one-loop integrals both analytically and numerically 

• routines for massive two-loop integrals 

• programs for numerical integration of two-loop diagrams. 

Here the graphical user interface and the text interface to Maple are presented. 

2.1 Introduction 

It is a well known fact, that high precision calculations of reactions between elementary 
particles cannot be done without the aid of computer programs, as the number of Feynman 
graphs can easily reach one hundred and exceed this as well. So there have been developed 
several different packages, each solving only one aspect or part of such calculations. More 
amazing is the fact that there is no package covering the whole procedure of calculating 
graphs. So someone starting such a high precision calculation is confronted with the fact, 
that he has to learn the syntax of several different programs. To overcome this problem, 
the development of XLOOPS, a program package covering all aspects of the calculation 
of Feynman graphs up to two-loop level was started. The following sections will discuss 
the structure of XLOOPS and introduce the graphical user interface (GUI), which makes 
XLOOPS an 'easy-to-handle' program package. 

^e-mail: bruecher@dipmza.physik.uni-mainz.de 
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2.2 The general structure of XLOOPS 



To avoid being tied to special computer systems or architectures, XLOOPS is designed to 
be as much portable as possible. Therefore 

. Maple Vi I, 

• TCL/Tk|]and 

• C++Q for numerical integrations 

were chosen as programming languages, as they are available for almost every computer 
architecture and all common operating systems, like Unix, Windows (3.1, 95 and NT) and 
VMS. 



File Model Mumber of Loops Extemail Legs 


Maple 








Helpi 












-8- 


-oo 





V/l IMbPLE V 

._|\| |/|_. Copyright (c) 1S81-1S90 by the University of Waterloo. 
\ tQPLE / All rights reserved, MAPLE is a registered trademark of 

< > Waterloo Maple Software 

I Type ? for help 

> LoopPath := Vusers/ul/iphuthep/TbruechecALDDps/MapleVRl/' : 

> read( Vusers/ul/iphuthep/lbruecher/SLoopa/ilapleVHl/loops. ma' ) : 
Warriifig: new definitiori for coKiJbine 

bytes U3ed=1014044, alloc=78e288, time=6 S33 

This is XLOOPS Version 1. Obeta test 

©1996 Lars BrTJcheij Joharjies Franskowski, Dirk Kreinier 



Figure 1: The main window of XLOOPS, which appears after starting the program 



This choice of languages can also be seen as an indicator for the program structure of 
XLOOPS. As already mentioned in the abstract XLOOPS consist of 

• a graphical user interface 

• routines for generating diagrams from particle input 
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• procedures for calculating one-loop integrals both analytically and numerically 

• routines for massive two-loop integrals 

• programs for numerical integration of two-loop diagrams. 

The graphical user interface should make the input of Feynman graphs as convenient as 
possible. This part is written in TCL/Tk and includes procedures for translating the user's 
input into commands suitable for the following parts. These parts, the second, third and 
fourth part, are written in Maple V. They provide procedures which insert Feynman rules 
and evaluate loop-integrals analytically as far as possible. The user might also use these 
procedures without the graphical front end. The fifth part is the numerical integration 
program for the remaining integrals in two-loop calculations which is automatically invoked 
if needed. This part will be covered by the talks of J. Franzkowski and A. Frink. 

When invoking XLOOPS the main window appears (see Fig. |^). It provides an overview 
of all topologies for a n-point m-loop function with n and m defined by the user in the main 
menu. By clicking on a topology the chosen graph appears in an extra window, suitable 
for inserting all virtual and non- virtual particles reacting (see Fig. Having inserted 
the particles one can choose whether the program should evaluate the result expressed in 
loop-integrals by clicking on the 'Evaluate' button or if the program should return also the 
loop-integral evaluated in logarithms and dilogarithm by clicking on the 'Evaluate Full' 
button. The program itself starts now the calculation of the graph by giving an appropriate 
command to the Maple V part of the package, which actually gives back the analytical 
result displayed in the Maple output section of the main window. If a two-loop calculation 
is performed and all numerical values are given, the program also starts the numerical 
integration of the remaining integrals. 

The result obtained by the evaluation can be saved or written out as a C program to 
a file. This is convenient for later use in a plotting program for example. For additional 
symbolic manipulations of the result, like inserting the renormalization conditions, a spe- 
cial window for direct Maple input exists. Moreover the program provides as additional 
features the possibility to insert the latest values from Particle Data Group Q for the 
particle properties as a menu entry. 

2.3 Examples calculated with XLOOPS 
2.3.1 The decay H+ — > W+h° 

To get more specific, the evaluation of — > W^h^ in the framework of the Two Higgs 
Doubletts Model 01 will be illustrated in this section. As an example we choose the top loop 
contribution to the self energy necessary for the wave function renormalization 
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and the triangle graph, which represents the actual one-loop correction 
discussion see 



(for a detailed 



mum 



File Model Mumber of Loops Extemail Legs Maple 

Standard Model i | 

♦ g Hlggs Douhletts Model | 



Evaluate Graph 




(q1«,0) 





Evaluate! Evaluate Full 




aose 







Evaluate Cr-aph 



(110,11) 



aid 

Help 



'A 




Evaluate 




Evaluate Full 




Close 



Figure 2: non-diagonal self energy graph and triangle graph contributing to W^h^ 

After choosing the '2 Higgs Doubletts Model' from the main menu and clicking on the 
graph we want to calculate, the window where the particles have to be inserted appears 
(see Fig. |2|). When finished with typing in the particles and clicking on 'Evaluate Full' the 
following result for the graph appears in the Maple output section: 



GO := [ 



ci = [1/32 



2 2 2 2 

e Mtop sin(alpha) cos (alpha) I (qlO - 4 Mtop ) 

s 

2 2 2 2 

sin(tw) Mw sin(beta) Pi 



2 2 2 2 2 2 2 

1/64 e Mtop sin(alpha) cos (alpha) (2 Ln(4 Pi MU ) I Pi (qlO - 4 Mtop ) - 2 



2 2 2 2 2 2 

Pi (-2 1 qlO +61 Mtop + Pi qlO - 2 Pi Mtop + I qlO gamma 
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2 2 2 2 2 

+ I qlO Ln(Pi) - 4 I Mtop gamma - 4 I Mtop Ln(Pi) - I Mtop Ln(Mtop ) 

2 2 

- I Mtop Ln(Mtop - I rho) 

2 2 

- I ("/.I (Ln(l - 7,1) - Ln(l + 7,1) + I Pi) - Ln(Mtop - I rho) - I Pi) qlO 

2 2 
+ 21 ("/.I (Ln(l - y.l) - Ln(l + 7.1) + I Pi) - Ln(Mtop - I rho) - I Pi) Mtop 

/ 2 2 2 4 

)) / (sin(tw) Mw sin (beta) Pi ) 
/ 

], CI] 

2 

- Mtop + I rho 

7.1 := Sqrtd + 4 ) 

2 

qlO 



After evaluating the tadpole contribution in a similar way one can insert the on-sheh 
condition, which in this case reads: 

^^hLo = ^2 ^^2 ^HOho{mlo) (9) 

This condition can be directly passed to Maple. In a similar manner the final result, after 
having summed up ah graphs, can also be exported as C language code. This is very 
helpful for producing plots as shown in |Q] . 

2.3.2 Flavour changing self-energy 

As a second example the flavour changing self energy Q from s ^ d in the framework of 
the Standard Model will be shown. 

After having inserted the particle properties from the main menu entry and having 
attached the appropriate particle at each particle line as shown in Fig. |3| the evaluation 
is again started with 'Evaluate Full'. In this case the program first starts the Maple part 
to evaluate the divergent part of the integral and the two-fold integral representation 
suitable for later numerical integration. As all numerical values are given, XLOOPS starts 
the remaining numerical integration. The result is displayed in the main window as a 
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function of the form factors. Each of the form factors' coefficients is displayed as a series 
in L» - 4. 




Figure 3: The self energy 



2.4 Availability and Outlook 

Currently a demo version of the program package, which does not include two-loop calcu- 
lations, is available at 

http : / /dipmza.physik.uni — mainz.de f^ ruecher / xloops.html 

where also additional information is accessible. The full version, which currently evaluates 
all one-loop graphs up to the 3-point function and all two-loop graphs up to the 2-point 
function, will also be accessible there, when it is fully tested. At the moment work to 
incorporate the 2-loop 3-point and the 1-loop 4-point function is in progress. 

In future, procedures that automatically draw all Feynman graphs for a given process 
and display them in Postscript format should be added as well as procedures for automatic 
renormalization. 
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3 



Automatic calculation of massive 
two-loop self-energies with XLOOPS 

Johannes FranzkowskiQ 



Within the program package XLOOPS it is possible to calculate self-energies up to the 
two-loop level for arbitrary massive particles. The program package - written in MAPLE 
||^, - is designed to deal with the full tensor structure of the occurring integrals. This 
means that applications are not restricted to those cases where the reduction to scalars 
via equivalence theorem is allowed. 

The algorithms handle two-loop integrals analytically if this is possible. For those 
topologies where no analytic result for the general mass case is available, the diagrams are 
reduced to integral representations which encounter at most a two-fold integration. These 
integral representations are numerically stable and can be performed easily using VEGAS 
ii- 

3.1 Structure of the XLOOPS package 

The aim of XLOOPS is to provide the user with a program package for complete evaluation 
of Feynman diagrams. The package consists of the following parts: 

• Input via Xwindows interface: 

In a window the user selects the topology which shall be calculated. A Feynman 
diagram pops up in which the particle names have to be inserted. 

• Processing with MAPLE: 

The selected diagram is evaluated. The necessary steps for reducing the numerator 
- the so-called characteristic polynomial - are performed by routines written for 
MAPLE. The result is expressed in terms of one- and two-loop integrals. 

''e-mail: franzkowski@dipmza.physik.uni-mainz.de 
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• Evaluation of one- loop integrals: 

One-loop one-, two- and three-point integrals are calculated analytically or numeri- 
cally to any tensor degree using MAPLE. 

• Evaluation of two-loop integrals: 

All two-loop two-point topologies including tensor integrals are supported by the 
MAPLE routines. For those topologies where no analytic result is known, XLOOPS 
creates either an analytic two-fold integral representation or integrates numerically 
with the help of VEGAS using C++. The other topologies can be calculated ana- 
lytically or numerically like the one- loop integrals. 

In the remainder of this report the MAPLE routines, especially the treatment of two-loop 
integrals is described in detail. 

3.2 The MAPLE part 
3.2.1 Input 

The user's input is inserted conveniently by using the Xwindows interface - there is another 
contribution by L. Briicher which describes this interface in more detail - but it can also 
be typed in directly in a MAPLE session if the conventions of the manual |p are respected. 

In any case the XLOOPS routines expect as input either a Feynman diagram or a 
single integral. In XLOOPS code Feynman diagrams look like 

EvalGraphl (3 , [bott om , charmbar , wp , bott ombar , gluon , 

gluon, charm, charmbar, bottom] ) ; 
EvalGr aph2 ( 6 , [up , upbar , gluon , upbar , up , gluon , 

gluon , up , upbar , up , upbar , gluon] ) ; 

EvalGraphl calculates one-loop, EvalGraph2 two-loop diagrams. The first argument de- 
notes the number of the topology, the second is the list of particles. Single integrals have 
the following notation: 

OneLoopSPt (pO)Pi5P2 >ql0,q20,q21 ,ml ,m2,m3) ; 
TwoLoop2Pt2(pO)Pii '''Oi ^1, s,qlO,ml ,m2,m3,m4,m5) ; 

The pi,ri,s determine the tensor structure Iq° l^^ l-^ and Iq" k^^ k^_l respectively - in 
the notation of section 2.3. The other arguments describe momentum components and 
masses. A detailed list of all conventions will be given in the XLOOPS manual the 
one- loop sector is also described separately 
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3.2.2 Characteristic polynomial 



In the case of a single integral XLOOPS skips the following steps. If a complete diagram 
has to be evaluated, the routines proceed as follows: 

• The topology - that means the information how the internal lines are connected - 
is generated. 

• The Feynman rules are inserted. At present XLOOPS knows all Feynman rules of 
the Standard model (electroweak and QCD) and its extension to two Higgs doublets. 
The incorporation of additional models is simple. 

• The symmetry factor of the diagram is determined. 

• The characteristic polynomial (the numerator of the diagram) is evaluated in terms 
of parallel and orthogonal space variables. For that purpose XLOOPS knows the 
rules for the SU(N) algebra and how to reduce strings of Dirac matrices. If necessary 
it takes the trace of the Dirac matrices. 

• Finally the code is expressed in terms of one- and two-loop integrals. 
3.2.3 Parallel and orthogonal space 

The evaluation of the Dirac algebra as well as of the integrals simplifies if one makes use 
of the splitting of the momentum components in parallel and orthogonal space variables 
- see D. Kreimer's contribution for details. The definition is simple: The parallel space 
describes the sub-space spanned by the external momenta, whereas the orthogonal space 
is the orthogonal complement of the parallel space. 

In the two-point case only one external momentum is present. Therefore the parallel 
space is one-dimensional. The loop momenta I and k are written as 




l-q 



(projection in direction of q) 



(the same decomposition holds for k) 



l-k 



loko — l±k± COS'!? 



The integration measure simplifies in the following sense: 
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one-loop case: 

oo 



D — 1 OO OO 

2Tr— 



dlo J dl± -2 

oo 



r(^) 

— uo 

The D-dimensional integral reduces to two one-dimensional integrations, 
two-loop case: 



oo oo 



J d^l jd^k= j dlo J dko j d^'H^. j d^-^k± 

oo 

O oo oo oo TT 

dlo J dko J dl± l^-'^ j dk± k^''^ j di} sin^^^ ^ 



— oo — oo 

D-1 D-2 OO OO oo 

47r 2 TT 2 



r(^)r(^) 

The 2Z)-dimensional integral is now replaced by five one-dimensional integrations. 



3.3 Two-loop integrals 

The aim of our two-loop routines is not only to solve special mass cases or kinematical 
regions. They also supply the general case where all internal masses are different and the 
external momenta can be completely arbitrary. Full tensor structure and not only the 
scalar case is supported. As a consequence these general routines cannot return analytic 
results for all topologies. Therefore we adopt the following strategy: 

1. separate the divergent parts (in any case analytically calculable) 

2. find a two- fold integral representation (in D = A) 

3. integrate numerically the two-fold integral 
In the following we comment these three items. 



3.3.1 UV divergent integrals 

For the separation of divergences it is necessary to find a convenient subtraction term. If 
a two-loop integral is multiplied by a term like the one in brackets 



/ 



/ 



Vi{l)V2{l)V3il + k)VA{k)V^{k) 



riii)r2ii)r3ii + k)mk)r5ik) V i\i + k) 
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it can be shown that the degree of divergence is decreased by m steps Q. To be specific, 
in the case of a logarithmic divergent integral which corresponds to the case m = 1 one 
gets a difference of two terms which turns out to be convergent. To recover the original 
integral it is necessary to add the subtraction term again. 

d^'ld^'k „ - / cPlfPk ^^^^ 



+ 



d^ld^k 



convergent 



i\i + kYVi{k)v^{k) 



divergent 



The first line which denotes the convergent part is treated in D = 4 and will be reduced 
to a two-fold integral representation which is solved numerically. The second line contains 
the divergent part which can be solved analytically m. D ^ A. 

3.3.2 Integration strategy 

We integrate directly in momentum space. With the splitting in parallel and orthogonal 
space variables one gets for two-point functions: 

oo oo TT 

j dH J d'^k = 8Tr'^ J dlodko Jlldl±kldk± JsinM^ 



analytically 

The Iq and ko integrations are left for numerical evaluation. All other integrations are 
performed analytically. 

Up to now we solved all two-point topologies for scalar and tensor evaluations. In the 
three-point case all topologies are solved for the scalar case. This will be subject of the 
contribution by A. Frink. At present our methods are applied to the three-point tensor 
and the four-point scalar integrals. 

In the following we concentrate on the two-point topologies: 
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The worst case is of course the master topology Q: 

ra jb uc ud 
'O ''± "'O ' 



dH / d^k 



l?il^,-^^kf.k'^^^ sin ^ 



CO oo oo oo 



= Svr^ J dlo J dko J dl^ J dk^ J d^ Vi---V5 

-oo -oo 

Vi = ll — — ml + ig 

V2 = ik + qf -ll-ml+ig 

^3 = [Iq + A;o)2 -lj_-k'{- 2lj_k_i_ cos i^-ml + ig 

Va = {ko - qf -k]_-m\ + ig 

'P5 = kQ — k'j_ — ml + ig 

The angular integration is elementary. For the orthogonal space integration the residue 
theorem is applied twice. The analytic result then is the following two-fold integral repre- 
sentation (written for the scalar case in dimensionless variables x,y). 



00 00 



dx J dy 



log(tt;3 + W2 + W5) + {similar terms} 
{wl - wi){wl - wl) 



Wl 



W2 



„ ml -ig 




X + \Y ^ — 



/, ,„ m\-ig 



i2 



WA = 



/ „ ml — ig 

3.3.3 Numerical evaluation 

The two-fold integral representation can be directly taken for numerical evaluation. We 
compared the results with available data. The data coincide perfectly with results from 
Shimizu, Kato and Fujimoto ||9|. They were also verified in several asymptotic limits |1C, 
0. 
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3.4 Future extensions 

The XLOOPS package is at present far from being saturated. The most important exten- 
sion will be the incorporation of two-loop three-point (and four-point) functions. In addi- 
tion a drawing routine for diagrams (generating Postscript output) will also be included 
as well as the possibility of evaluating complete processes and the option to renormalize 
the processes automatically. 
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4 



Massive two-loop vertex functions 

Alexander FrinkQ 



Calculating massive two-loop vertex functions by splitting integrations into parallel 
and orthogonal space components has been demonstrated to give a convenient twofold 
integral representation suitable for numerical evaluation. This method is now extended 
to more topologies, including the crossed topology, graphs with divergent subloops and 
infrared divergent diagrams. The connection between this representation and physical and 
anomalous thresholds is examined. 

4.1 The crossed vertex function 

The scalar crossed vertex function (Fig. ^) can be calculated similarly to the planar vertex 
function (Fig. |^) with the orthogonal/parallel space integration technique presented in 
||^]. The volume element is written as d^/ d'^k = ^dlo dko dli dki ds dt da dzl\/\ — z^, where 
^0, ^1, /co and k\ are the components of the loop momenta parallel to the external momenta, 
s = /j^, t = fcj^ and z is the cosine of the angle between /_|_ and k^. The procedure can be 
divided into the following main steps: 




Figure 4: Topologies covered in this article 

• linearization of propagators in l\ and k\ by shifting — *■ + ^ij ^ + k\ 

• integration over orthogonal space angles (a and z) 
^email: Alexander. Frink@Uni-Mamz.DE 
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• integration over li and ki with Cauchy's theorem, avoiding cuts resulting from the 
z integration 

• integration over s and t with Euler's change of variables 

• numerical evaluation of the remaining integrations over Iq and in a finite region 

The main complication arises here as both loop momenta / and k have to flow through two 
propagators in common. To perform the z integration as in we have to apply a partial 
fraction decomposition in these two propagators, in order that in the difference term the z 
dependence drops out. Each term of the partial fraction decomposition can be calculated 
separately and similarly to the planar vertex function as in Q, however the condition in 
which way to close the contour for the li and ki residue integrations is different for both 
terms: Iq + ko — ei^O and /o + + £2 ^0 resp., where ei and 62 are components of the 
external momenta. As a consequence poles from the difference term which lie on the real 
axis also have to be taken into account. These poles give rise to several unbounded areas 
in the {Iq, k^) plane in addition to triangles. Nevertheless these add up to zero outside the 
finite region shown in Fig. ^. 




Figure 5: Effective integration region in the (Iq, ko) plane for the crossed vertex function 
After the residue integrations we have the intermediate result 

1 n 00 00 

10 . . ^ 



Vc = YlJJ^^o dko jds Jdt Cj 



A , 

1 1 



dijt + bij + cijC d2jt + b2j + C2jC 

(10) 



a3jt + b3j + C3jC (^ajt + bj + CjsY - 4st 



where Cj, a^j, b^j, c^j, aj, bj and cj are rational functions of and ko, and Aj is a subset 
of the area in Fig. ||. Some of the dkj and Ckj vanish. Euler's change of variables for s 
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and t reduces the problem to integrals of the forms / ln(2;^ + rx + s) /{x^ A-px + q) dx and 
/ arctan(x + a)/{x'^ + px + q)dx, which can be expressed in terms of dilogarithms and 
Clausen functions. Further details can be found in 

As for the planar vertex function a correlation between the various coefficients in 
Eq. (10) and physical thresholds resulting in an imaginary part of the diagram can be 
found: if an dkj or Ckj is zero, then the accompanying bkj corresponds to a two-particle 
cut. The coefficients bj correspond either to three-particle cuts or the two-particle cut 
q'^ > (mi -|- 771-2)^ which is somewhat hidden due to the partial fraction decomposition. 



4.2 Anomalous Thresholds 



The connection between the coefficients in our two-loop three-point representation and 
physical thresholds has been pointed out in for the planar and above for the crossed 
topology. Now it will be demonstrated that anomalous thresholds are also accounted for 
correctly in this representation. As an example, let us consider the one-loop three-point 
function in Fig. |^. By searching for solutions of the Landau equations HI, one can find a 
singularity which corresponds to no physical threshold at 1 -|- 2fiifj.2H3 — //^ — /^l ~ Ms = 
where fij = (mf + 1712 + rn^ — mj — p'j)mj/ (2mim2m3). If we arbitrarily choose all masses 
equal, mi = 1712 = = 1, and pi = 6, P2 = 5, these equations are fulfilled for P3 = 
—4 — \/l5 ~ —7.873. A plot of the real and imaginary part of this diagram for different 
values of in the vicinity of this point is shown in Fig. ^(left). 

Let us now turn to the two-loop example Fig. |^. The Landau equations give us a 
singularity for this graph at the same values as above for arbitrary 1114, and mg. We 
recall the four- fold integral representation for the planar graph 

vp = y ffdiodkoc f - — ^ — TT f — — ^ = (11) 

^JJ ^ ^ J (t + to)it + t'o) J s + so J (at + b + csY - Ast ^ ' 

J-^ Tj 

to and t'o are quadratic functions in k^. As was pointed out in real roots of to and 
tg inside the integration triangle correspond to physical thresholds p\ > (m2 + m^)'^ and 
P2 > {m^ + TTT-i)^ respectively. These conditions are fulfilled here. Let us denote the roots 
of to and I'q as and k^ resp. Then it can be shown that below the anomalous 

threshold k^^^ < k^^^' and above k^^^ > k^^' . However, for the numerical integration over 
Iq and /cq (after partial fraction decomposition) the program only has to make a distinction 
between t^^ > (real part only) and tg^ < (additional imaginary part). Merely numerical 
stability gets slightly worse in the direct vicinity of the threshold. Fig. ^(right) shows the 
expected plot of real and imaginary parts. 



25 



real part 
imaginary part + 



< >000000000000t>00000«» 



+++++++ 

Y 



real part 
imaginary part 



<j@oooooooooooooooooo< > 



-7.883 -7.878 -7.873 -7.868 -7.863 -7.883 -7.878 -7.873 -7.868 -7.863 



Pi 



Pi 



Figure 6: Plot of real and imaginary parts of a vertex function in the vicinity of an 
anomalous threshold, one-loop (left) and two-loop (right) 

4.3 Divergent diagrams 

The method for calculating massive two-loop vertex functions by splitting integrations into 
parallel and orthogonal space components works without modifications only for convergent 
integrals. The complication stems from the non-elementary integration over the cosine of 
the angle between l± and k±, z, in D = A — 2e dimensions 



1 



z"^ A + Bz + i7] 



dz 



Z^) 2 



1 



A + Bz + ir] 



dz . 



(12) 



which can only be expressed in terms of hypergeometric functions instead of a simple 
square root. 

In order to avoid this, one has to subtract a simpler function with the same structure of 
divergence. The difference is finite and can be calculated with the parallel/orthogonal space 
technique in four dimensions, whereas the subtracted function should be calculable in D 
dimensions. A first example for an UV divergent scalar vertex function at zero momentum 
transfer has been given in |^] . This will now be extended to other genuine vertex topologies 
as well as infrared divergent diagrams. 



4.3.1 Infrared divergences 

As an example let us consider the planar graph depicted in Fig. |^ with qf = ml, = "^i 
and niQ — > 0. We choose the loop momentum I to flow through Pi . . . P3 and k through 
P4 . . . Pq. This graph exhibits an infrared divergence for /c — > 0. By subtracting the same 
diagram, but with P3 replaced with P3 = PsI^^Q' power-counting for small k is improved 
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by one unit, so the difference is now finite and can be calculated in four dimensions: 



This involves the same main steps as for the crossed vertex function (see above). Since 
infrared divergences are endpoint singularities, all steps can be applied individually to both 
parts of the integrand, and the infrared divergence will not show until the last numerical 
integration over feg near /cq — > 0. But there the subtraction guarantees a well-defined limit. 

Let us now demonstrate the integration steps for the subtracted term. The z integration 
is trivial, since none of the propagators depends on z. Therefore there is no cut in the 
integrand which has to be avoided during the residue integrations. We have the freedom 
to close the contours analogously to the original term i.e. for /q + /cq > in the upper 
and for /q + ^0 < in the lower half plane. Then the integrations over li and ki give 
conditions of the form Iq^x and ko^y for the individual propagators to contribute, which 
results in triangular regions listed in Table ||. The overlap of these triangles coincides with 
the overlap of the triangles from 



Table 1: Non- vanishing triangles from the subtracted term 



propagators conditions 



(^1= 


P^) 


k 


+ ko<0 


lo + ei- qz 


> 


ko + e2 + Qz 


> 


(Pi, 


Pe) 


lo 


+ ko<0 


lo + ei- Qz 


> 


ko>0 




iP2, 


Pa) 


lo 


+ ko>0 


/O - 62 - Qz 


< 


ko-ei + qz 


< 




Pe) 


lo 


+ ko>0 


^0 - 62 - Qz 


< 


ko<0 




(Ps, 


Pa) 


lo 


+ ko>0 


lo<0 




ko-ei + 


< 


iP3, 


P5) 


lo 


+ ko<0 


lo>0 




kQ + e2 + qz 


> 



The analytic evaluation of the subtracted part involves calculating one-loop three- 
point functions up to 0{£) This subtraction procedure is also applicable for small but 
non-zero mg in order to improve numerical stability and to extract the leading logarithm. 

4.3.2 Graphs with divergent subloops 

Ultraviolet divergent scalar diagrams as in Fig. ^ can also be calculated by subtracting an 
appropriate quantity with the same structure of divergence. This can be done by setting 
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the masses in the inner subloop to zero. Let k flow through the divergent subloop and I 
in the outer loop. Defining P4 = P4|^^^q and P5 = -Psl^g^o' have 



Performing the z, li and ki integrations as above, we obtain the intermediate result 

2 



Vuvfimte = ^ JJdlo dko j ds JdtC 



(s + so)(s + s'q) 



J=i 

1 1 



Viat + b + cs)'^ -Ast J^at + b + csY-^st 



(15) 



This can be readily expressed as a twofold integral representation over dilogarithms. The Zq 
and kQ integrations extend again over triangles. Note that the boundaries of the triangles 
depend solely on the external momenta and not on the internal masses. 

The calculation of the subtracted part is done by first integrating the massless subloop 
in k, giving a result proportional to (^^)^^. The evaluation of the one-loop three-point 
function with non-integer powers of propagators can be done at least numerically. 

The topology in Fig. |^ can be reduced to the topology in Fig. |^ by partial fraction 
decomposition in the propagators P3 and Pg- Since the integration region in the {lo,ko) 
plane is equal for both terms, they can be added again before integrating over Iq and ko. 
The case ms = mg is handled by differentiating Vjjv with respect to m^. 



4.4 Outlook 

Current work is done in the direct calculation of tensor integrals. A general method using 
subtraction terms similar to above has been proposed by D. Kreimer Q. 
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